Optical Phonon Anomaly in Bilayer Graphene with Ultrahigh Carrier Densities 
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Electron-phonon coupling (EPC) in bilayer graphene (BLG) at different doping levels is studied 
by first-principles calculations. The phonons considered are long-wavelength high-energy symmetric 
(S) and antisymmetric (AS) optical modes. Both are shown to have distinct EPC-induced phonon 
linewidths and frequency shifts as a function of the Fermi level Ef- We find that the AS mode has 
a strong coupling with the lowest two conduction bands when the Fermi level Ep is nearly 0.5 eV 
above the neutrality point, giving rise to a giant linewidth (more than 100 cm -1 ) and a significant 
frequency softening (~ 60 cm" 1 ). Our ab initio calculations show that the origin of the dramatic 
change arises from the unusual band structure in BLG. The results highlight the band structure 
effects on the EPC in BLG in the high carrier density regime. 

PACS numbers: 63.20.kd; 63.22.Rc; 73.22.Pr; 78.30.Na 



Electron-phonon coupling (EPC) is an important effect 
in monolayer graphene (MLG) pj. Interesting phenom- 
ena such as the renormalization of the phonon energy 
Q, the Kohn anomalies Q, and the breakdown of the 
adiabatic (Born-Oppenheimer) approximation 0,01 have 
been reported. It is expected that even more intriguing 
effects should be found in AB-stacked bilayer graphene 
(BLG) where both the band structure and the doping 
level (i.e., the Fermi level Ep) can be tuned through 
the applied electrical gates [fMI, allowing for the con- 
trol of a delicate interplay between electrons, phonons, 
and photons flol - flS] ]. Previous investigations mainly fo- 
cused on situations with charge carriers near the charge 
neutrality Dirac point E D (i.e., \E F - Ed\ <0.4 eV) fiol — 
12 . IToT ]. In this energy range, phenomena such as the 
phonon mode renormalization, Raman broadening, and 
the phonon frequency shift induced by EPC have been 
successfully predicted by a simplified tig ht-binding (TB) 
model and effective mass theory 

Recent progress in fabricating electrolytic and 
ionic-liquid gates provides the possibilities of doping 
MLG and BLG with ultrahigh charge-carrier densities of 
\n\ > 10 14 cm" 2 and of tuning the Fermi level close to 
the van Hove singularity (VHS) point at M in the first 
Brillouin zone (BZ) [lilHij- Such a high carrier density 
will embark interesting technological applications includ- 
ing supercapacitors [2l| . transparent electrodes 1221 . and 
high performance organic thin film transistors [23(. In 
this regard, understanding the carrier dynamics will be 
a crucial step for the potential electronic device applica- 
tions. Distinct many-body effects and superconducting 
instability have been observed in doped graphene when 
the Fermi energy approaches the VHS point [20j . Since 
superconductivity also occurs in graphene-related sys- 
tems such as graphite-intercalation compounds (GICs), 
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FIG. 1: (Color online) Band structures of bilayer graphene 
calculated from (a) a simplified tight-binding model with 
only nearest-neighbor (NN) coupling and (b) density func- 
tional theory. The x axis is in units of 2n/ao, where ao is 
the graphene lattice constant. All possible transitions (I- VI) 
through the phonon modes are labeled in (b). The neutrality 
point Ed has been shifted to energy zero. 



the study of EPC in BLG might shed new light on the 
underlying mechanism of superconductivity in GICs (23 - 

Raman and infrared (IR) spectra are powerful tools to 
probe the EPC in MLG and BLG by providing the infor- 
mation of the phonon linewidth and frequency shifts as 
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a function of doping and field strengths \4 
Although the aforementioned models [13, [14( are known 
to work well for low-energy charge carriers, deviations 
between theoretical predictions and experimental obser- 
vations become significant when \Ep — Ep>\ > 0.4 eV [271 ]. 
Such a discrepancy calls for a full consideration of band 
structure effects from first principles, as the electronic 
structure in this regime is expected to be considerably 
modified from its low energy part. 
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In this work, we show that the band structure of BLG 
at high doping levels plays a critical role in the EPC of the 
two long-wavelength high-energy symmetric (S) and an- 
tisymmetric (AS) optical modes in the high carrier den- 
sity regime. For both modes, the linewidth 7 exhibits 
a dramatic dependence on the electron and hole doping 
level. In particular, the AS (E u ) mode exhibits a giant 
EPC-induced phonon linewidth (more than 100 cm -1 ) 
and significant softening (60 cm -1 ) when the Fermi level 
Ep is 0.5 eV above the neutrality point En. Despite 



an intensive study of the EPC in BLG [13|, \l£ l28( , this 
feature has not been reported so far. This giant enhance- 
ment originates from the strong coupling between the AS 
mode and two conduction bands. The strong EPC in this 
high carrier-density regime might affect the performance 
of BLG-based devices. 

The first-principles calculations reported in this work 
are performed using the Quantum ESPRESSO code 
29l 30| . The electronic structure is obtained using the 
local density approximation (LDA) within the density- 



functional theory (DFT), and the core- valence inter- 
action is modeled by norm-conserving pseudopotentials 
[1H . The wave functions of the valence electrons are ex- 
panded in plane waves with a kinetic energy cutoff of 70 
Ry. A vacuum region of 20 A has been introduced to 
eliminate the artificial interaction between neighboring 
supercells along the z direction. The relaxed C-C bond 
length is 1.42 A and the interlayer distance is 3.32 A for 
BLG. 

The phonon frequencies and associated eigenvectors 
were computed using the density-functional perturbation 
theory (DFPT) [i^l , details of which have been presented 
in our previous work [H, The self-energy n q „(u;) of 
a phonon with wave vector q, branch index is, and fre- 
quency u q „ provides information on the renormalization 
and damping of that phonon due to the interaction with 
other elementary excitations. Following the Migdal ap- 
proximation, the self-energy induced by the EPC in BLG 
reads O]: 



n q „H = 2 
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where e TO k is the energy of an electronic state Imk) with 
crystal momentum k and band index m, /(e m k) the cor- 
responding Fermi occupation, and r\ a positive infinites- 
imal. For a given mode uj = loq, the phonon linewidth 
is 7 = — 2Im(n q i / (wo)) and the phonon frequency shift is 
Aw = ^[Re(n qi/ (cj )|.E F -n q „(iJo)|.E J ,=o]- 

The EPC matrix element in Eq. (1) is given by 



3m n ( k , q) 



(2) 



where SV sc f = l/s C /(w q ) — V sc f(0) is the variation of 
the self-consistent potential field due to the perturba- 
tion of a phonon with wave vector q and branch index 
v. Variations of the potential field SV sc f are calculated 
through self-consistent calculations to find the poten- 
tial field for both perturbed and unperturbed systems. 
The perturbed phonon mode is handled with the frozen- 
phonon approach [33j. The DFT calculations have been 
carried out on a dense 201x201 fc-grid within a minizone 
(0.4x0.4) enclosing the BZ corner K in the reciprocal 
space. This is equivalent to 500x500 fc-grid sampling in 
the whole Brillouin zone. Finally, the EPC matrix el- 
ements are computed using Eq. (2). By changing the 
Fermi level Ep in Eq. (1), we can investigate the depen- 
dence of 7 and Alj on different doping levels, assuming 
the EPC matrix elements are unchanged. This approxi- 
mation is justified by the small dependence of the EPC 



matrix elements on doping for the T phonon modes in 
graphene 34|. For all the linewidths calculated below, 
we used a parameter 77 = 5 meV. 

In Fig. [T]we show the band structure of BLG obtained 
by the TB model (with only nearest-neighbor (NN) inter- 
actions) and DFT calculations. For BLG, the low-energy 
dispersions within a TB model can be well described by 
E(k) = ±t 2 /2 ± V^/4 + ^|/(fc)| 2 , with t x = -2.7 eV 
and ^2 = 0.36 eV, the intralayer and interlayer hopping 
parameters, respectively [35] . The low-energy parabolic 
dispersions predicted by the TB model are in agreement 
with the DFT result. Nevertheless, there is still a notice- 
able difference: In Fig.[ljb), the band dispersions have an 
evident asymmetry between the conduction and valence 
bands, while the TB model shows nearly symmetric band 
structure relative to the neutrality point [Fig.[TJa)]. This 
electron-hole asymmetry for the low-energy charge car- 
riers has already been revealed by infrared spectroscopy 
0] and becomes more significant for the high-energy car- 
riers. In particular, our first-principles calculations pre- 
dict that there is a crossing between the two conduction 
bands (cl, c2) at about 0.1x27r/ao from K along the K- 
M direction in the BZ [3q |. The energy is around 1.0 eV 
above the neutrality point, as shown in Fig. [Tfb) . This 
feature is not captured by the simple TB model, indicat- 
ing the necessity of first-principles calculations. As we 
will show later, the band structure plays a crucial role in 
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FIG. 2: (Color online) Calculated linewidths 7 for the (a) 
symmetric (S) and (b) antisymmetric (AS) mode as a function 
of Ef- Note that the vertical scale has been changed in (b). 
The neutrality point Ed is set to be energy zero. Insets are 
schematic plots of the S and AS modes. 
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FIG. 3: Calculated frequency shift Auj for the (a) symmet- 
ric and (b) antisymmetric mode as a function of Ef. The 
neutrality point Ed is set to be energy zero. 



understanding the physics of high-energy charge carriers. 

The asymmetric band structure of BLG has impor- 
tant implications on the phonon linewidth and frequency 
shift because the S and AS modes couple with different 
bands in BLG [33]. As indicated in Fig. [IJb), the al- 
lowed transitions [i.e., nonzero EPC matrix elements in 
Eq. (2)] for the S mode are transitions I (vl-cl), II (v2- 
c2), V (vl-v2), and VI (cl-c2). In contrast, the AS mode 
only couples with transitions III (vl-c2), IV (v2-cl), V 
(vl-v2), and VI (cl-c2). When E F is tuned to be 0.5 
eV above the neutrality point, the transition VI becomes 
active because the energy requirement is satisfied. On 
the other hand, the crossing of the two conduction bands 
(cl, c2) in this energy range provides a larger electronic 
phase space in the Brillouin zone for the phonon mode 
to couple with, as will be shown later. As a result, the 
linewidth of the AS mode is expected to be significantly 
enhanced. 

In Figs. Uta) and [U[b), we show the calculated 
linewidths for the S and AS modes with respect to 
Ep. The S and AS modes (Hujq ~0.2 eV) considered 
in this work have been schematically shown in the insets. 
The corresponding frequency shifts Aw are presented in 
Fig. GHa) and[«2[b), respectively. 

In the low doping regime with \Ep — Ep>\ < hwo/2 ~ 
0.1 eV, the S mode can be in resonant coupling with 



the electron-hole pair from the the top valence band 
(vl) and the bottom conduction band (cl). As a re- 
sult, the linewidth of the S mode is constant, as depicted 
in Fig. [2j^a). In contrast, the AS mode has a negligi- 
ble linewidth in this range. Our calculations of phonon 
linewidth and frequency shift in this energy region rc- 
jroduce the features found in previous DFT calculations 



28 1 and agree reasonably well with the experimental data 
111 ]. This is also consistent with a previous study based 
on a continuum model [l3j |. 

In comparison, dramatic differences are found for the 
ultrahigh electron doping when Ep — Ep> > 0.5 eV. The 
calculated 7 for the AS mode is significantly larger than 
that of the S mode. For example, when Ep — Ep, = 
0.70 eV, 7 is 150 cm -1 for the AS mode, around 20 
times larger than that of the S mode (8.0 cm -1 ). Corre- 
spondingly, there is a significant softening (~ 60 cm -1 ) 
for the frequency shift when Ep —0.6 eV, as shown in 
Fig-E^b). Further increase of the Fermi level will result in 
an even more significant increase of the phonon linewidth 
(not shown). In contrast, the phonon softening becomes 
smeared out as Ep increases. Although such a doping 
level requires large electron densities, it may be achieved 
either by directly applying electrical gate in experiment 
13] or by chemical doping [U . For Ep — Ep, = 0.5 eV, 



the desired electron density in BLG is estimated to be 
n ~ 6xl0 13 cm -2 , still an order of magnitude smaller 
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FIG. 4: Calculated joint density of states (JDOS) N(E f ,lj ) 
as a function of the Fermi level. The neutrality point Ed is 
set to be energy zero. 



than the doping limit achieved in monolayer graphene 
(~ 4xl0 14 cm -2 ) [16j. Therefore, we expect that this 
high electron doping regime will be realized in experi- 
ment |17j and that the giant phonon linewidth as well as 
the significant phonon mode softening may be verified by 
Raman or infrared measurements. 

This giant linewidth enhancement mainly arises from 
the large phase space associated with the high doping 
range. In Fig.[4]we plot the joint density of states (JDOS) 
N(E f ,oj ) = & Ew\fM-fM]5[eiv-e kj .+fi Uo ] 
as a function of Ep. The JDOS increases dramatically 
when the Fermi level is in the range of Ep — Ed >0.5 
eV, which allows more electronic states to couple with 
the phonon modes. The profile of N(Ep) is similar to the 
phonon linewidth profile of the AS mode, indicating that 
the giant enhancement is mainly an electronic structure 
effect. 

In summary, our first-principles calculations find that 
the phonon linewidths and frequency shifts for the 
long-wavelength high-energy optical modes in bilayer 
graphene exhibit a distinct dependence on the electron 
and hole doping due the intriguing interplay between the 
unique band structure and the phonon modes in this sys- 
tem. In particular, we predict that the linewidth for 
the antisymmetric mode could be significantly enhanced 
when the Fermi level is tuned to be 0.5 eV above the 
neutrality point. 
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